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Abstract 

We have calculated the Bjorken-x dependence of the kaon and pion valence 
quark distributions in a statistical model. Each meson is described by a 
Fock state expansion in terms of quarks, antiquarks and gluons. Although 
Drell-Yan experiments have measured the pion valence quark distributions 
directly, the kaon valence quark distributions have only been deduced from 
the measurement of the ratio Uk{x) /u^x). We show that, using no free 
parameters, our model predicts the decrease of this ratio with increasing x. 
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1. Introduction 

A determination of parton distribution functions (PDFs) in mesons is 
important as a test of QCD. A Drell-Yan experiment has measured the pion 
valence quark distributions [1] , but the kaon valence quark distributions have 
only been deduced from the measurement of the ratio uk(x) /u n (x) [2]. The- 
oretical calculations of pion PDFs have used the Dyson Schwinger equations 
(DSE), the Nambu-Jona-Lasinio (NJL) model, instantons, constituent quark 
models, and statistical models. For a recent review see Holt and Roberts [3] . 
Kaon PDFs have been calculated with the NJL model [U E], a meson cloud 
model [6], a valon model [7j. and the Dyson Schwinger equations [H[9]. 
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In section 2 we summarize the statistical model, review our calculation 
for pion PDFs, and describe our calculation for kaon PDFs. In section 3 
we compare our valence quark ratio to experiment and other theoretical 
calculations. 

2. statistical model 

Zhang and collaborators [TUl HTJ H2] have used a simple statistical model 
to calculate parton distribution functions in the proton. They considered 
the proton to be an ensemble of quark-gluon Fock states, and used detailed 
balance between states to determine the state probabilities. A Monte-Carlo 
program was used to generate momentum distribution functions for each 
state, from which proton PDFs were determined. Their model, using no free 
parameters, predicted an integrated asymmetry (d — u) = 0.124, consistent 
with the experimental value 0.118 ±0.012 measured in deep inelastic scatter- 
ing (DIS) [13]. Their results were also in good agreement with the d(x) — u(x) 
distributions measured in DIS [H] and Drell-Yan experiments [T5| ITS]. 

2.1. statistical model for the pion 

We have used Zhang et al.'s statistical model to calculate PDFs for the 
pion that are in good agreement with experiment and other calculations [T7] . 
We assumed a light sea of uu, dd and gluons. The Fock state expansion for 
the pion is 



with i the number of uu pairs, j the number of dd pairs and k the number 
of gluons. The leading term in the expansion, {ijk} = {000}, represents 
the valence quark state ud. The probability of finding the pion in the state 
\{ud}{ijk} > is fHjk = \cijk\ 2 , with 




(1) 



i,j,k 




(2) 



Detailed balance between any two Fock states requires that 



PijkN(\{uud}{ijk} >— > \{uud}{i'j'k'} >) 
Pi'j'k'N(\{uud}{i'j'k'} >->■ \{uud}{ijk}) >, 



(3) 
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in which N(A — > B)) is the transfer rate of state A into state B. Transfer 
rates between states are assumed to be proportional to the number of partons 
that can split or recombine. Taking into account the processes q -H- qg, 
g «-> q q, and g gg, 

Pijk _ 1 rr 2 + 2i + 2j + n 

Pooo ~ i\{i + l)\]\{3 + l)\l} Q (2 + 2i + 2j){n + l) + r ^ ) 

This equation, together with the normalization condition (|2]), determines the 
Pijk- The 7r + sea is flavor symmetric, since the equation is symmetric in i 
and j. 

The average number of partons in the pion, n n , is given by 



[ [u(x) + d(x) +u(x) + d(x) +g(x)]dx = 4.74 . (5) 
Jo 



For the pion we assumed massless partons, and used the Monte Carlo 
event generator RAMBO [18] to determine the momentum distribution f n (x) 
for each n-parton state, n = 2 + 2{i + j) + k. The flavor distributions for 
each Fock state are 

dijk{x) = fn(x)(l+j) , Uijk(x) = fn(x)(l + i) , (6) 

Uijk{x) = fn{x)l , dijk(x) = fn{x)j , (7) 

and for the gluons 

9ijk(x) = fn{x)k . (8) 

The PDFs are found by summing these distribution functions over all values 
of {ijk} 

u(x) = ^2p ijk u ijk (x) = d(x) , (9) 

i,j,k 

d{x) = ^Pijkdijkix) = u(x) , (10) 



and 



g( x ) = ^2 Pijk 9ijk(x) . (11) 
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The valence quark distribution function is 

v(x) = u(x) — u{x) = d(x) — d(x) (12) 

and the momentum distribution functions satisfy 
1 

x[u(x) + d(x) + u(x) + d(x) + g(x)]dx = 1 . (13) 

In order to compare our valence quark distributions to experiment, we 
carried out an evolution in Q 2 . We determined the starting scale of our 
distributions by requiring that the first and second moments of our valence 
quark distribution at Q 2 = 4 GeV 2 be equal to those found by Sutton et al. 
[T9"] . This gave us a starting scale of Ql = 1.96 GeV 2 . We used Miyama and 
Kumano's code BF1 [20] for the evolution. Our results, shown in Fig. [TJ are 
in reasonable agreement with experiment. The non-zero values at x ~ 1 are 
contributed by the first term in the Fock state expansion, which consists of 
two massless partons, for which the momentum distribution is constant. 
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Figure 1: Our results (solid curve) for the valence quark distribution xv(x) in the pion, 
compared to the experimental results of Conway et al. PQ. The dashed curve shows our 
results without any evolution, which correspond to a scale of Qq = 1.96 GeV 2 . The solid 
curve shows our results evolved to Q 2 — 16 GeV 2 of the E615 experiment. 
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2.2. statistical model for the kaon 

As for the pion, we assume that the kaon has a light sea of uu, dd and 
gluons. The Fock state expansion for the K + is 

\K + >=J2d jk \{us}{ijk}>. (14) 

Again, assuming that transfer rates between states are assumed to be pro- 
portional to the number of partons that can split or recombine, and including 
the processes q -H- q g, g -H- q q, and g -B- gg, 

Pijk = 1 t-t 2 + 2i + 2j + n 

pooo i\{i + l)\]\ 2 AA( 2 + 2z + 2j)(n + l) + ^' 

We find the average number of partons in the kaon Hk = 5.23. The K + sea 
is asymmetric, with < uu >= 0.424 and < dd >= 0.685. As in the case of 
the proton, the valence u quark provides more pathways for annihilation of 
the u quarks in the sea, whereas the d quarks can only annihilate on other 
sea quarks. Detailed balance then requires an excess of dd over uu. 

For the kaon we also used RAMBO to determine momentum distributions 
for each n-parton state. If all partons are considered to be massless, including 
the s, PDFS for the kaon differ from those of the pion, because uk > n^. The 
kaon's momentum is shared among more partons than the pion's momentum, 
so the momentum fraction contributed by its valence quarks 

f xu K {x)dx = 0.46 (16) 
Jo 

is less than the momentum fraction contributed by the pion's valence quarks 

xu n (x)dx = 0.50 . (17) 



o 



This is seen in the valence PDFs shown in Fig [2} The momentum distribu- 
tions at x ~ 1 are determined by the probability of finding the meson in the 
leading term of its Fock state expansion, the valence state, which is 0.14 for 
the pion and 0.10 for the kaon. 

We also considered the case in which the s was given its current quark 
mass M = 100 MeV, but all other partons were considered massless. We 
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Figure 2: Valence quark distributions xv(x) for the kaon (solid curve) and the pion (dashed 
curve), assuming massless partons. The PDFS are calculated at our starting scale of 
Ql = 1.96 GeV 2 . 



then have two distributions, f n M(%), for the massive s, and f n o(x) for all the 
other partons, considered massless. We determined PDFs for the kaon, and 
compare them to the pion, in Fig. [3j As expected, the uk{x) distribution 
shifts to lower x, and the sk{x) distribution peaks at higher x. Parton 
numbers and momentum fractions for the pion and the kaon are shown in 
Table [B 

Table 1: Parton numbers and momentum fractions for the pion and the kaon, calculated 
at the starting scale of = 1.96 GeV 2 . For the kaon, < x > is the momentum fraction 
calculated assuming a massless s, and < x >m is the momentum fraction calculated for 
an s mass M = 100 MeV. 
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Figure 3: Valence quark distributions xu 7r (x) for the pion (solid curve), xuk(x) (short- 
dashed curve) and xsk{x) (long-dashed curve) for the kaon, with an s mass M = 100 
MeV. The PDFS are calculated at our starting scale of Q\ = 1.96 GeV 2 . 



3. Comparison with experimental ratio and other theoretical cal- 
culations 

Badier et al. [2] determined the valence quark ratio uk{x)/u^{x) from 
Drell-Yan experiments with K~ and 7r _ beams incident on a platinum target. 
By symmetry this ratio is equal to the valence quark ratio uk(x)/u^(x) 
for the positively charged mesons we have discussed above. We compare 
our results to experiment in Fig. |4j Evolution has very little effect on 
the ratios. The massless parton calculation agrees better with experiment 
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Figure 4: The valence quark distribution ratio uk{x)/u v (x). The left panel shows our 
calculation for massless partons. In the right panel the calculation used M — 100 MeV 
for the s. The dashed curves were calculated for the starting scale; the solid curves are 
for distributions evolved to Q 2 = 27 GeV 2 for comparison with the Badier et al. [5] 
experiment. 
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than the calculation that included the s mass. However, as seen in Table [TJ 
the momentum fraction carried by the Uk(x) distribution, drops from 23% 
(massless s) to 20% (massive s) and the momentum distribution shifts to 
lower x, as seen in Fig. |3j Both effects cause the ratio to increase for x < 0.2 
and decrease for larger values of x. In Fig. [5] we compare our results to other 
calculations, and to experiment. 
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Figure 5: Calculations of the valence quark distribution ratio u K +(x)/u 7r +(x) compared 
to experiment [2]. The solid line is our calculation using M = 100 MeV for the s mass. 
The long dashed curve is our calculation using massless partons. The short-dashed line 
is the DSE prediction of Nguyen et al. [8]. The dot-dash curve is the NJL calculation of 
Davidson and Arriola [5] The dotted curve is the valon model calculation of Arash [7] for 
a massive s. 

The best agreement with experiment is the recent DSE calculation of 
Nguyen et al. [8], which used a full Bethe-Salpeter amplitude. Other calcu- 
lations predict an x-dependence similar to ours. Davidson and Arriola, using 
an NJL model, found the same trend of a decreasing ratio with increasing 
x, as did the earlier NJL calculation of Shigetani et al. J4]. Arash, using a 
valon model with equal masses for the u and s valence quarks in the kaon, 
found reasonable agreement with experiment, but for a massive s, the ratio 
did not agree as well. These models do not fit the experimental values of 
the ratio for the mid-range of x as well as the DSE calculation. The models 
predict a wide range of values for the ratio as x — > 0. Data for the ratio in 
the region x < 0.2 is needed to test the models. In their review article, Holt 
and Roberts [3] have noted that the experimental data for the ratio are 'not 
of sufficient quality to test and verify our understanding of pion and kaon 
structure'. It is important to make higher statistics measurements of both 
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the pion and kaon parton distributions, and to extend the measurements to 
lower x. 

4. Conclusions 

We have used a simple statistical model, developed for the calculation of 
parton distribution functions in the proton, to calculate the parton distribu- 
tion functions of the pion and the kaon. We find that the ratio of valence 
quark distributions, uk(x)/u w (x), shows the expected decrease with increas- 
ing x. 
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